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Abstract 


By using an integral formula, we easily calculated the von Neumann 
entropy of any quantum state. 


1 Introduction 


The expression for von Neumann entropy|1]-[2] is 
S(p) = —TrpInp. (1) 


where p is the density matrix of any quantum state. Now, we introduce an 
integral formula 


1° 2 
If IneIn(1+ $ )dz = alna- a, a>0. (2) 
T Jo ax 


Since any density matrix p is semi positive definite, without loss of generality, 


we get 
2 


1% p 
—plnp = —p — — ln zln(1 + —)dz. 
plnp p = nz In( + 7a) de (3) 
Tracing both sides of Eq. (3), we have 
S(p) = —Trplnp (4) 


1 o0 p 
= -Trp- F InaTr In(1 + = )dz. 
T JO 


x 
According to the formula Tr ln A = Indet A, Eq. (4) can be changed as 
2 


1 co 
S(p) = -Trp — If ln zlndet(1 + a )dz, (5) 
TJo 


r2 
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which is exactly the integral representation of von Neumann entropy S(p), where 
det (1 + 2) is a Fredholm determinant in fact. For a Fredholm determinant of 


n multiplied by n matrix, there is the following expansion formula: 
I+tTrA +t detA, n=2 
I+tTrA+Ẹ (TrA)? — TrA?] +t det A, n=3 


det(I+ŁA) = I+tTrA+Ẹ4[(TrA)? — TrA?] 
+Ë (TrA)? — 3(TrA)(TrA?) + 2(TrA3)] +t4detA, n=4 
31 
(6) 
Then, for n = 2, 
Trp? det p? 
det(1 + 5) = a ~ (7) 
= a Trp? = 
z i 72 
For n = 3, 
2 T 2 T 2 2 _T 4 det 2 
dt = rp, |(Trø) ro]  detp (8) 
x x Qa4 gô 
-j Tro [(Trp*)? — Trp*] det p\” 
= 2 a z3 
For n = 4, 
p 
det(1 + 72) (9) 
pa e , are - Tro] 
7 5 ge! Q!a4 
_[(Prp?)? — 38(Tr@) (Trø) + 2Prp%)] , (det p\* 
3!a8 x4 ` 
It is already known that any pure state satisfies p? = p, so 
Trp = Trp? = Trp” = 1. (10) 


According to Eq. (5), we show that for n by n density matrix p for a pure state, 


IOS 2 
S(p) = -Trp— 1y ln zlndet(1 + Pda (11) 
T JO T 
17s 1 (det p\” 
= -1-2/ mami+ = + ( : e) jdz. 
T Jo x x” 
The result of Eq. (11) is very concise and beautiful. 
For any pure state, because 
det p? = (det p)? = det p, (12) 
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so, 
det p = 0,1. (13) 


As we all know, det p can only be 0 for any pure state, thus its von Neumann 
entropy S(p) = 0. 
2 Eigenspectrum expression of S(p) 


A similar transformation P can always be found to diagonize the density matrix 
p. Let the diagonalized matrix be D, then 


PtP = D = diag(\1, A2,---,An); (14) 
where A; is the eigenvalue of p. We have 
p? p? 
det(1 +=) = det[PT (1 + —)P] (15) 
x x 
D? 
= det(1+ =) 
= [[a+ me 
k=1 
then 
1° p 
S(p) = —Trp-—— ln zlndet(1 + —)dx (16) 
T Jo x 
1 hae = Az 
= -l- =i tarte Oa 2 da 


1% ° A? 
-1--S°/ Inzin(l+ S)dz 
Set 0 xv 


—1— S$ (Ag ln àk — Ax) 


k=1 
-X Ayn An, 
k=1 


which is consistent with the traditional expression of S(p). 


3 Series expansion of S(p) 


Usually, the computation of solving the eigenvalue of density matrix p is very 
large, so we try to find a simpler way to calculate von Neumann entropy S(p) 
to reduce the computation. Based on the form given by Eq. (4), we speculate 
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whether it is possible to rewrite S(p) into some series form. Firstly, we should 
consider 


n(1 + se = 5! a , (17) 


then 


= (18) 


2 oo kp p2k 
_ (-1)*Trp 
Trln(1 4 72) J 


However, unfortunately, by substituting Eq. (18) into Eq. (4), we found that 
the integration is divergent. To make up for this deficiency, we try to write the 
integral in Eq. (2) form in piecewise form and use the substitution method to 
perform some deformation treatment on the integral, which means 


o0 2 
f lngln(1 + <3) de (19) 
1 2 co 2 
= i Ina In(1 + Z da +f ln gzln(1 + de 
0 T 1 x 
i ‘na In(1 + a?r? 
= i Inv In(1 «dx | ania = 2 PA 
0 x 0 x 


However, the first integral of Eq. (19) after written as the series form is still 
divergent, so we have to continue to transform it into the following form 


1 2 
J IneIn(l + $ )dz (20) 
0 x 


1 
f In x In(x? + a°)dz — 4 
0 


1 z2 
= | lngln(1 + ~ )dx +2lna -— 4. 
0 a 
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Thus, 


oo 2 
a 
Ing In(1 + qa de (21) 


i 
1 2 1 27,2 
= I lnzln(1 4 ave f Mp tal aa y + 2Ina-—4 
0 a 


x2 


1 oO x? 1 
—%)' lnag =ar? l 
= -f T i a+ fo 3 xr+2lna— 4 


k=1 
© (_g2)k pl _a—2)k pl 
= 5i “) | Ing « g? pet lng: zdr +2Ina—4 
k=1 0 k=1 o 
oo (—a7)* 3 (=a-")" 
= X +2Ina—4 
zI 2 
= k(2k — 1) E k(2k +1) 
© (_]1)k 2k —2k 
= > ) [ = 5 2 z| +2lna— 4 
— k (2k — 1) (2k + 1) 
We notice that a~?* = ya —a?*)! for a < 1 and that series C ai 
1=0 k=1 
and Da T -i Gkr are all convergent for a < 1, so Eq. (21) can be changed 
into T form as follows 
ee) a2 
i Inv In(1 + ade (22) 
0 
= ak 1 S 2kyl 
1 + 21 4. 
- >< i TE arip ao eee 
k=1 1=0 
Now, let a be p, then 
S(p) 1 =f Trini +2) (23) 
= -l--— n n = 
p A xTr qa de 


1% 1)" 1 = 2kyl Trp 
= Tr(1 
2 k rP i= ap- 
2Indetp 4Tr1 
Gy l T 


In this way, we simplify the calculation of von Neumann entropy S(p) to the 
problem of solving the series of Trp?* and Tr(1 — p?*)!, which gives us the 
possibility to calculate von Neumann entropy S(p) by using numerical methods. 
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4 Discussion 


In fact, for the problem of integral divergence caused by the series expansion of 
2 
In(1 + $+) in equation (2), we can also use Borel resummation method to solve 


it. Let 
(ert on 


S = ln(1 + x) = 3 = Xar, (24) 


-(—x)*t+! 


zyr By using 


then az = 


S 


aan ant 2 
cont X /_ k+1 

f e 5 xt) dt 

0 ! 


II II 
= 
pee 
eos 
+ 8 
Ghose Di 
~ 
ai 
a. 
a 
~~ 
a 
~ 


and applying this method to Eq. (2), we have 


Lf? a? 
-/ Inen(. + Sade (26) 
TT 0 xv 


1 foe) (s ome ar a ere 

= -+f ma f £ i dtdx 
T Jo i Ca (k +1) 
1% et a2t 

= -+f ma f f (e — 1)dtdx 
T JO 0 
1 co. -t [oe) 

= -1f =| Ina(e” =? — 1)dxdt 
TJo t Jo 
1 fet [> a2t 


II 

| 

l 
= 
5 
a 
Fie? 
Pi 
Y 


— 1)dzdt. 


However, this solution is not conducive to solving the problem of von Neumann 
entropy, but rather brings considerable complexity to the problem. 
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